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Abstract 

An analytic approach is developed to predict the shape and displacement with voltage 
in the quasi-static limit of round LaRC Thunder Actuators. The problem is treated with 
classical lamination theory and Von Karman non-linear analysis. In the case of classical 
lamination theory exact analytic solutions are found. It is shown that classical lamination 
theory is insufficient to describe the physical situation for large actuators but is sufficient 
for very small actuators. Numerical results are presented for the non-linear analysis and 
compared with experimental measurements. Snap-through behavior, bifurcation, and 
stability are presented and discussed. 

Introduction 

Thunder actuators are devices constructed from an isotropic laminate of aluminum, 
LaRC SI adhesive, piezo-electric PZT, and a metal backing of either steel or brass. The 
LaRC adhesive is a 1 mil thick solid thermoplastic and the aluminum layer is indented to 
provide electrical conductivity after the actuator is cured at 250 deg C. After the curing 
process they form a bow shape after reaching room temperature. If simply supported or 
arranged in a clamshell arrangement they produce a linear motion when an external 
voltage is applied. 

A simple method for controlling the manufacturing process and predicting the 
actuation amplitude has long been sought for Thunder actuators. This has been recently 
been treated by Campbell[3,4] in his treatment of rectangular actuators. However, the 
round actuators are particularly important in the manufacture of pumps and speakers. 
Because of the non-linear nature of the problem, the solution is not necessarily a dome. 
Unlike the case with the rectangular actuators, this is a highly undesirable effect. What is 
needed is a method for predicting the shape before the device is manufactured. 

One approach to this problem is to use classical lamination theory. This is a very old 
subject based on the theory of composite laminates. With laminates, temperature effects 
are very important. As a result a very rich theory has been developed to minimize these 
temperature effects. Although the theory is sufficient to find an optimal lay-up to 


1 



minimize temperature effects, it is quite deficient in predicting the actual shape in all 
cases. The reason for this is because it is a linear approximation to a problem that is 
essentially non-linear. As a result, a non-linear correction to the theory is needed to 
describe large actuators. 

The type of correction used here is one that is valid for large displacement and small 
strain. This correction is usually attributed to Von Karman [2] in plate theory. 
Campbell[3,4] has applied this to rectangular actuators with good success. This work was 
inspired by Hyre[5] who used this technique in his treatment of the temperature effects of 
orthotropic laminates. He used this in combination with a simple Raleigh-Ritz energy 
method to obtain solutions exhibiting the snap-through behavior. Isotropic laminates 
produce domes that approach cylinders when the laminate is large. It is possible to treat 
the round actuator case in exactly the same way as the rectangular case by converting to 
cylindrical coordinates before calculating the energy. 

Solution by Classical Lamination Theory 
Derivation of equations 

Classical lamination theory [1] assumes strains are small, and amplitudes are shallow. 
In addition, it assumes both shearing and extensional strains perpendicular to the surface 
are zero. These assumptions follow the Kirchhoff-Love hypothesis for shells. 

Under deformation the thin plate is described by three variables - u, v and w, where u 
is the displacement in x direction, v in the y and w is the height of the plate in the z 
direction. The variables u, v and w are themselves functions of x, y and z. It is assumed u 
and v follow the relation. 


u = u Q - z 

V = v 0 - z 


d_% 
3 x ’ 
3 w 0 
3 y ’ 


( 1 ) 


where the 0 subscript denote the quantity is measured about the middle surface. The 
strains are given by [2], 
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The results of Equation 1 and Equation 2 imply, 
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where k x , Ky, and K X y are the x, y and twist curvatures. The middle surface strains and 
curvatures are given by. 
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The relationship between stress and strain is given by, 
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where Q is the reduced stiffness and Q is the transformed reduced stiffness. For isotropic 
media. 
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where Y is Young's modulus and v is Poisson's ratio. To find the force and moments we 
first find the stresses for the k'th layer by, 
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Now we use the relationship, 
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where h is the thickness. By Equation 9 and 8 we find, 
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where A, B, and D are the extensional, coupling and bending stiffness matrices and are 
given by, 


(id 

D «= 3 l.feWS -*-')• 

where the z's are measured with respect to z=0 (at a point half way through the 
thickness). Since each layer is isotropic the symmetry of the stiffness matrices will be the 
same as the reduced stiffness matrix. Moreover, since in this case the layers are not 
symmetric about z=0 the coupling stiffness will be non-zero so that, 
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We now have everything we need to find the thermal loading and actuation loading so 
we may solve the problem. 

Thermal Loading 

The thermal loading takes two forms. The first is extensional and the second is 
bending. The problem is greatly simplified since we have only isotropic layers with equal 
thermal strains in each direction. The thermal loading may be calculated by, 
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Actuation Loading 

By analogy with thermal loading, 
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where t is the thickness of the PZT, V is the applied voltage, d 3 i is the actuation constant 
in a plane perpendicular to the poling axis and, 


N A =d 3 fa+Q^V, 
M A =d 3l (Q l + Ql(z 3 +$V. 


(16) 


Here it is assumed we only have actuation on the third layer. 

Thermal and Actuation Displacement 

To find the resulting strains and curvatures in terms of the actuation loading we first 
add both the thermal and actuation sources as, 
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Now we use the results of Equation 10 so that, 
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Inversion of the above gives, 



By using the results of 12 and inverting the 6x6 matrix we find, 
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(19) 


( 20 ) 


Y -iy ~ K atv = 0 ' 

So we have thus found the curvature in terms of known quantities. The shape is that 
of a dome. The simplicity of the relationship comes about because the layers are 
isotropic. The solution is valid for small actuators such as the round ones used for pumps 
and speakers. 

Solution by the Von Karman Non-Linear Approximation 

The Von Karman approximation assumes large displacement and small strain. This 
suggests including second order terms in w in addition to the linear terms for u, v and w. 
Expanding Equation 2 to include these terms give, 
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This combined with the results of Equation 1 give, 
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The total energy takes the form, 

U= j J j (dU 0 -dW T -dW^, ( 23 ) 

where dU 1 is the stored elastic energy volume density and dWj and dWA are the 
temperature and actuation contributions. The elastic energy portion of it for the k'th layer 

is. 
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For the actuation portion we have, 


(dW^ = (o A x e x + a$ e) 3 , 

= (<2i + G 2 ) 3 ^31 y B .v+ (2i + 2 2 ) 3 ^3i y e / 


(26) 


Now that we have the general form for the energy integral we may try a Raleigh-Ritz 
type solution. The method involves making a guess for the solution. This guess should 
look qualitatively like the expected solution yet be flexible enough to allow for 
adjustment. This adjustment comes in the form of changing the guess based on a finite 
number of parameters which the guess contains. The parameters are adjusted in such a 
way as to minimize the total energy of the system. In this situation it would be 
advantageous to choose a guess that is close enough to the classical lamination result that 
in the linear limit (small scaling) the solution approaches the classical lamination result. 
Such a guess takes the form, 

w° = j (a x 2 + b y 2 ). (27) 


Classical lamination theory assumes no shearing strain between layers for thermal 
expansion. If we make the same assumptions here we find u*-* and v^ must take the form, 

u° = c x- ^ a 2 * 3 - a b x y 2 , 

(28) 

v° = d y-X a 2 f - \ a b x 2 y. 
o ‘ 4 

The ansatz represented by Equation 27 and 28 has been used successfully by 
Campbell[3,4] in his treatment rectangular actuators and Hyer [5] in his treatment of 
orthotropic laminates. Equations 22, 27 and 28 give, 

s x =c-X a bf-az, 


8 v = d - ^ a b x 2 - b z. 


(29) 


Y.<y = 0 - 

All that is needed now is to use Equation 23 with Equations 24, 25 and 26 and 29 to 
find the total energy. To do the round actuators, however, we must first convert to 
cylindrical coordinates. The transformation is, 
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( 30 ) 


x = r cos (0 + cp) 
y = r sin (0 + cp) 
dx dy dz = rdrdQ dz 

The angle cp is an arbitrary phase angle that comes about because of symmetry 
arguments. The total energy is obtained by integrating r from 0 to R, and 0 from 0 to 2jt . 
z is integrated through the thickness. The solution is obtained by minimizing the total 
energy with respect to a, b, c and d. The result of the minimization is, 

D l a + K A b - B x c - B 2 d + K l b 2 + 2 K l a b - K 2 b c - K 2 b d + 

K 3 ab 2 + M'= 0 

K 4 a + D l b-B 2 c-B l d + K l a 2 + 2K l cib-K 2 ac-K 2 ad + 

K 3 a 2 b + M' = 0 (31) 


- B ] a - B 2 b + A l c+A 2 d-K 2 ab-N' = 0 

- B 2 a - B l b + A 2 c + Aj d - K 2 a b - AT = 0, 

where, 

K, = ^{B l + B)^ 

k 2 =±(a 1+ a^r 2 

(32) 

K * = m ( 3 A > +A ^ 

K 4 =D 2 + f# 


The A's, B's and D's have the same meaning as before and N' and M' is the same as 
well. The radius R is the length and width of the actuator. It is a simple matter to show 
that the solution in the linear limit (small scaling - R=0) matches the classical lamination 
result if we take a=-K x , b=-Ky, c=s x and d=Ey. The non-linear equations are coupled and 
third order so it is unlikely that an exact analytic solution exists. It may be possible to 
find close approximate solutions to them, however. In any event they may be solved 
numerically with no difficulty. 
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Unlike the classical lamination result the non-linear analysis predicts unequal 
curvatures that depend on the magnitude of the scaling and the aspect ratio. It also 
predicts multiple possible solutions. The nature of these solutions are dome-like solutions 
which approach cylinders in the limit of large scaling. 

As a comparison for a particular case we compute a simple example. For this case we 
choose a 3 layer actuator with 1 mil LaRC Si, 6.8 mil PZT and 2 mil brass with R=.6 in. 
The results are shown in Figures 1 and 2. Although we have two solutions, because there 
is also an arbitrary phase angle, 4>, these are actually the same solution rotated by 90 
degrees. In fact, any rotation of one of those solutions is also a solution. If R is small 
enough we have a single dome solution. The solutions shown in Figures 1 and 2 are 
sometimes referred to as a saddle mode but technically they are not true saddles since the 
curvatures in each direction have the same sign. As we shall see there are some cases 
where we can actually have 3 stable solutions. This is different from the perfectly square 
actuator where there are at most 2 stable solutions. 



Figure 1. Non-linear Solution #1. 



Figure 2. Non-linear Solution #2. 
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Figure 3. Non-linear Solution #3. 


Stability 

The stability of each solution is determined by whether or not the solutions are 
maxima or minima. If a particular solution represents a minimum then it is stable. If it is 
a maximum, it is unstable. The total energy of the system may be viewed as a function of 
several variables. As such, once an extremum is found, the test whether it represents 
minuma is determined by whether the determinant of F is positive definite where, 

Det[F] >0, 

(32) 

F - d 2 U 
1 j d Xj d x j 

Fortunately d U/d xj is represented by the left hand side of Equation 31 (to within a 
multiplicative factor) so all we have to do are the xj differentiations to find the matrix 
elements. 

Bifurcation 

As stated earlier, The solution may take on different shapes depending on the scaling. 
To demonstrate this, numerical results were completed for a square piece having the same 
material properties as before except we vary the size from R=0 to R=2 in. For a small 
piece the shape is a dome as predicted by classical lamination theory. At some critical 
size the shape turns suddenly into a dome with unequal curvatures. In the limit of large 
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scaling it becomes a cylinder. Figure 3 is a plot of both the x and y curvatures. In this 
case the critical size is at approximately 1 inch square. 



Figure 3. Bifurcation plot of x and y curvature as a function 
of size for a round actuator. 


This result shows that there are 3 possible solutions past the critical size to thickness 
ratio. Experiment shows at a large enough size the perfect doming is rarely seen in 
Thunder actuators. We feel this is due to slipping during cure - that there is a slight 
preferential direction in the strain. In Rainbow actuators, on the other hand, the doming 
solution is seen at all sizes. In the case of Rainbow slipping is not possible because the 
base layer is created by a chemical reducing process. 

Perfect doming solution 


A single equation may be derived for the perfect doming case. This is obtained by 
solving Equation 31 in the simple case where a=b and d=c. The result is, 
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Expressed in terms of the dome height this becomes, 


(34) 
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This has an exact analytic solution. It is, however, quite complicated and will not be 
presented here. 

Deflection under non-inertial distributed load 


To estimate the deflection under load we write the total energy as, 


U = 
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\dU n -dW T -dW 
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2 Jt r w P dr 
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where w is the height as a function of r and P is the pressure. Following the same procure 
as before we find, 
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In terms of the dome height, 
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(37) 


(38) 


This result impliies a higher stiffness than actually occurs. By taking the isotropic and 
linear limit and comparing with exact plate theory it can be shown this result is stiffer by 
about 4/3. The reason for this is because the system is heavily over constrained. The 
second order polynomial used here can't accurately account for the change in shape under 
load. For approximation purposes, however, it is sufficient. 
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